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Abstract — In this paper, we consider the downlink of large- 
scale multi-cellular OFDMA-based networks and study perfor- 
mance bounds of the system as a function of the number of users 
K, the number of base-stations B, and the number of resource- 
blocks TV. Here, a resource block is a collection of subcarriers 
such that all such collections, that are disjoint have associated 
independently fading channels. We derive novel upper and lower 
bounds on the sum-utility for a general spatial geometry of base 
stations, a truncated path loss model, and a variety of fading 
models (Rayleigh, Nakagami-m, Weibull, and LogNormal). We 
also establish the associated scaling laws and show that, in the 
special case of fixed number of resource blocks, a grid-based 
network of base stations, and Rayleigh-fading channels, the sum 
information capacity of the system scales as 6(73 log log K/B) 
for extended networks, and as 0(731oglog K) and f2 (log log if) 
for dense networks. Interpreting these results, we develop some 
design principles for the service providers along with some 
guidelines for the regulators in order to achieve provisioning 
of various QoS guarantees for the end users and, at the same 
time, maximize revenue for the service providers. 

I. Introduction 

With the widespread usage of smart phones and the increas- 
ing demand for numerous mobile applications, wireless cellu- 
lar networks have grown significantly in size and complexity. 
Consequently, the decisions regarding the deployment of base 
stations, the ratio of the number of subscribers to the number 
of base stations, the amount to be spent on purchasing more 
bandwidth, and the revenue model to choose have become 
much more complicated for service providers. Understanding 
the performance limits of large multicellular networks and 
the optimal balance between the number of base stations, the 
number of subscribers, and the bandwidth to achieve those 
limits are critical components of the decisions made. Given 
that the most significant fraction of the performance growth 
of wireless networks in the last few decades is associated [1] 
with the cell sizes and the amount of available bandwidth, the 
aforementioned issues become more important. 

To that end, in this paper we analyze the achievable down- 
link information rate in large multicellular OFDMA systems 
as a function of the number, K, of users, the number, 73, 
of base-stations, and the number, TV, of available resource- 
blocks. Here, a resource block is a collection of subcarriers 
such that all such collections, that are disjoint have associated 
independently fading channels. The contributions of this paper 



can be summarized as follows. 

• For a general spatial geometry of the base-stations and the 
end users, we develop novel upper and lower bounds on the 
average achievable rate as a function of K, 73, and TV. 

• Then, we consider two asymptotic scenarios in network size: 
dense networks and extended networks in which (user) nodes 
have a uniform spatial distribution. We evaluate our bounds 
for Rayleigh, Nakagami-m, Weibull, and LogNormal fading 
models along with a truncated path-loss model. To evaluate 
the bounds, we utilize various results from the extreme value 
theory. We also specify the associated scaling laws in all 
parameters. 

• With the developed bounds we consider four different 
scenarios. In the first scenario, we consider afemtocell network 
and develop an asymptotic condition for K, 73, and TV to 
guarantee a non-diminishing rate for each user. In the second 
and third scenarios, we consider extended multicell networks 
and derive bounds for the choice of K/B, i.e., the ratio of 
the number of users to the number of base stations, in order 
for the service provider to maximize the revenue per base 
station and at the same time keep the per-user rate above a 
certain limit. We analyze two different revenue models: for the 
overall service the users are charged (1) on the number of bits 
they are served; (2) a constant amount. Finally, we consider an 
extended multicell network and develop asymptotic conditions 
for K, B, and TV to guarantee a minimum return on investment 
for the service provider. 

Calculation of achievable performance of wireless networks 
has been a challenging, and yet an extremely popular problem 
in the literature. The performance of large networks have been 
mainly analyzed in the asymptotic regimes and the results have 
been in the form of scaling laws [2]— [ 11] following the seminal 
work by Gupta and Kumar [2]. Unlike these studies, our main 
bounds are not asymptotic and we take into account both a 
distance based power-attenuation law and fading into account 
in our model. Scaling laws for channel models incorporating, 
both, distance based power-attenuation and fading have been 
considered in [12], [13]. We assumed a truncated path-loss 
model unlike these works, which assume an unbounded path- 
loss model. Unbounded path loss models affect the asymptotic 
behavior of the achievable rates significantly. For instance, 
the capacity scaling law of 6(log7l) found in [12] arises by 
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exploiting infinite channel-gain of the users close to the base- 
station, whereas, without path-loss the scaling law changes to 
B^oglogif). Further, our analyses take take into account the 
bandwidth and number of base-stations in large networks. 

II. System Model 

We consider a multi-cellular time-slotted OFDMA-based 
downlink network with B base-stations (BS) and K active 
users as shown in Fig. 1. The base-stations lie in a disc of 
radius p — R (p > R > 0), and their locations are arbitrary 
and deterministic. The users lie in disc of radius p and their 
locations are uniformly distributed within this disk. We denote 
the coordinates of BS i (1 < i < E) with (a,i,bi), and the 
coordinates of user k (1 < k < K) with ix k ,y k ). 




Fig. 1: OFDMA downlink system with K users and B base- 
stations. 

We now describe the channel model. We assume that the 
OFDMA subchannels are grouped into resource blocks, across 
which the BSs schedule users for downlink transmission. We 
denote the complex-valued channel gain over resource-block 
n (1 < n < N) between user k and BS i by /i" fc . We assume 
that h™ k is composed of the following factors: 

K k = mzz<k, (i) 

where /3R~ k denotes the path-loss attenuation, 

Ri^k = max{r , \J {x k - a.;) 2 + (y k - h) 2 }, (2) 

/?, a (> 1), ro (< R) are constants, and the fading factor 
v™ k is a complex-valued random variable that is i.i.d. across 
i,k,n. Therefore, for any given user-BS combination, the 
channel-gain over each resource-block is independent 1 of the 
channel-gains over other resource-blocks [14]. Currently, we 
keep the distribution of uf k general and specific scenarios 
will be studied in subsequent sections. Assuming unit-variance 
AWGN, the channel-SNR between user k and BS i across 
resource-block n can now be defined as 

llk = Kk\ 2 = (3 2 R7 2 k a W?J 2 - (3) 

'Note that this assumption is for simplicity and our analyses can be easily 
generalized to the case with dependencies across resource-blocks. 



We assume the base-stations have perfect coordination 
among each other via a backhaul network and they allocate 
resource-blocks and downlink powers in each time-slot jointly 
such that the achievable sum-rate of the system is maxi- 
mized. We denote the corresponding (sum-rate maximizing) 
user scheduled by BS i across resource-block n by Uf and 
the corresponding (sum-rate maximizing) power allocated by 
P". We assume that, in each time-slot, every BS allocates 
powers to scheduled users subject to a sum-power constraint 
P con . Hence, the set of feasible user allocations and power 
allocations are: 

U = {< V i, n : 1 < < < K V i, n}, and 

V 4 KVi,n:£y<Pcon V (4) 

n 

Assuming the availability of perfect channel-state informa- 
tion of all users' channel-gains at every BS, the achievable 
sum-rate of the system can be written as 

B N , _n pn v 

= EE^+ 1+y ^V J . (5) 

i=l n=l ^ 1 + 7 i.^ ^ ' 

= max log ( H — 1 , (6) 

where x := {x k V k}, y := {y k V k}, u := {v™ k V i, k,n}, 
U := {Up V i,n}, and P := {P™ V i,n}. 

Note that C x . y ^ v (U ,P) is a random variable, which is a 
function of the BS-locations (aj,bj) for all i (deterministic), 
and the random variables (xk,y k ) and v™ k for all (i,k,n). 
In the following section, we derive novel upper and lower 
bounds on the mean of C x . y ^(U, P) to determine the scaling 
laws and develop our network-design guidelines. To state the 
scaling laws, we use the following notations. For two non- 
negative functions f{t) and git): 

1) f(t) = Oig{t)) means that there exists a positive con- 
stant c\ and an real number v\ such that jit) < Cigit) 
for all t > ri. 

2) f{t) = f2(p(t)) means that there exists a positive con- 
stant C2 and an real number r-i such that fit) > czg(t) 
for all t > r 2 . In other words, g(t) = 0(f{t)). 

3) f{t) = e(s(t)) means that f{t) = 0(g(t)) and 
fit) = Qigit)). Note that if lim^oo = C3 for some 
constant C3, then fit) = 0(g(i)). However, vice versa 
is not always true. 

III. Bounds on Achieved Sum-Rate 

In this section, we derive performance bounds for the system 
model defined in Section II. The expected achievable sum-rate 
of the system, using (5), can be written as 

C*=E{C a , V j,(U,P)}, (7) 

where the expectation is effectively over the SNRs 

{7&,Vi,*,n}. 
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The following theorem uses extreme-value theory [15] to 
derive novel upper and lower bounds on (7). 

Theorem 1. The expected achievable sum-rate of the system, 
C* , can be bounded as follows: 



log (1 + P con max fc 7"^ 



< C* < min j E { log (1 + Peon max 7^)}, 



N 



(8) 



Proof: Proof is given in Appendix A. ■ 

IV. Application of Bounds to Networks 

The bounds in Theorem 1 are quite general and can be 
applied to a variety of network and channel settings. In the 
following subsections, we consider the specific cases of dense 
and regular extended networks, and apply (8) to obtain the 
performance bounds and associated scaling laws. 

A. Dense Networks 

Dense networks contain a large number of base-stations that 
are distributed over a fixed area. Typically, such networks 
occur in dense-urban environments and in dense femtocell 
deployments. Thus, in our system-model, dense network cor- 
responds to the case in which p is fixed and K and B grow. 
The following theorem builds on Theorem 1 to find bounds 
on achievable sum-rate of the system for Rayleigh fading 
channels. 

Theorem 2. For dense networks with Rayleigh fading down- 
link channels, i.e., v™ k ~ CA/"(0, 1), 

( log(l + P con l K ) + 0(l))BNf™(r, B, N) 

< C* < (log(l + PcopZa-) + 0(l))BN, (9) 

where r > is a constant, Ik = /3 2 rQ 2a log —j^f-, 

fu (r, B. N) = t--tz — ^ — r „„ -i n , rrrr, and a, a are 

the mean and standard deviation of \v™ k \ 2 (p = a = 1 for 
Rayleigh fading channels). The following scaling laws result 
from (9): 

C* = O(BNloglagK), and 
C* = n(min{B,N} log logiC). 



(10) 



Proof: The complete detailed proof is provided in [16]. A 
detailed sketch with three intermediate lemmas (without proof) 
are given in Appendix B. To summarize, the first lemma, i.e, 
Lemma 1, uses Cantelli's inequality and Theorem 1 to show 
that 

C* > /r(r,P,7V)^E{log(l + P con max 7 ^)}. 

The second lemma, i.e, Lemma 2, finds the distribution 
of channel-SNR 7™ fc under Rayleigh-distributed \vfy.\ and a 



truncated path-loss model defined in (1). The third lemma, i.e, 
Lemma 3, uses Lemma 2 and extreme-value theory to show 
that (maxfc 7" fe — Ik) converges in distribution to a limiting 
random variable with a Gumbel type cdf, that is given by 



cxp( 



7/3 2 



), X € (—00, 00). 



(11) 



Thereafter, we use Lemma 1, Lemma 3, and [9, Theorem A. 2] 
in (38)-(49) to obtain the final result. 



Corollary 1. If we have lo s^ N ^> ^ then a tighter upper 
bound for Rayleigh-fading channels can be stated as C* = 
0(BAriogl2£^V). 

The proof sketch is given at the end of Appendix B. For 



log KN 
N 



» 1 



detailed proof, see [16]. Note that the condition 
is easy to satisfy since, typically, the number of users scale 
much faster than the number of resource blocks (particularly 
when the number of resource blocks are fixed). 

Similar results under different fading models are given by 
the following theorem. 

Theorem 3. If \v™ k \ belongs to either Nakagami-m, Weibull, 
or LogNormal family of distributions, then, for dense 
networks, the scaling laws for the upper bounds are 

For Nakagami-(m, w): C* = 0(BN log log K) 

For Weibull(X, t): C* = 0(BN log log^ K) 
For LogNormal '(a, J): C* = 0(BNy/\ogK). 

and the scaling laws for the lower bounds are 

ForNakagami-(m,w).- C* = Cl(BNf^ N (r, B, N) log log if) 

For Weibull (A, t): C* = n(BNf° N (r, B, N) log log? K) 
For LogNormal (a, u): C* = n(BNf£ N (r, B, N)^/TogK). 

Proof: Proof sketch given in Appendix C. For detailed 
proof, see [16]. ■ 



Corollary 2. // '-^^ » 1, ^f^- » 1, and » 1 

for Nakagami-m, Weibull, and LogNormal distributions, re- 
spectively, then tighter upper bounds on the scaling laws can 
be found as follows: 

For Nakagami-(m, w): C* = O (BN log 1 -^^ L ^j 



log t KN 



For Weibull (A, t): C* = Ol BN log ^ 



For LogNormal (a, uj): C* = Oi BN log 



, logT KN 
N 

v/log K N ' 



N 



The proof sketch is given at the end of Appendix C. For 
detailed proof, see [16]. 

B. Regular Extended Networks 

In extended networks, the area of the network grows with 
the number of nodes, keeping node density constant. Here we 
study regular extended networks, in which the base-stations 
lie on a regular grid. We assume the grids are hexagonal as 
illustrated in Fig. 2 and the distance between two neighboring 
base-stations is 2R. Hence, the radius of the network p = 
Q{RVB) for large B. 
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Base Stations 



Fig. 2: A regular extended network setup. 



The following theorem gives the performance bounds, and 
the associated scaling laws for regular extended networks and 
Rayleigh fading channels. 

Theorem 4. For regular extended networks with Rayleigh 
fading downlink channels, i.e., iy" k ~ CAf(0, 1), 



>g(l + Pcon^A') + 0(l))BNf™(r,N) 
<C* < (log{l+P C0n l K ) + O(l))BN, 



(12) 



where l K = pr^loggfe, f^{r,N) 



(l+r 2 )~V 



CO 



P ft 2 T 2 - 2a 

iconl-> < o 



V3(2q-2) 

channels), . The associated scaling laws are: 



standard deviation of \v™ k \ (p, 



N+(/j.+r<T)c ' 

and p, a are the mean and 
(7 = 1 for Rayleigh fading 



C* = 0(BN log log ^) 



c* = n siog log 



K 
B 



(13) 



Proof: Proof sketch is provided in Appendix B. For 
detailed proof, see [16]. ■ 

Corollary 3. If ° N B 3> 1, then a tighter upper bound on 

/ log KN . 

the scaling law can be written as C* = 0[BN log — 

This condition is easy to satisfy since, typically, the number of 

users scale much faster than the number of resource blocks. 

The proof sketch is given at the end of Appendix B. For 
detailed proof, see [16]. Note that for fixed number of resource 
blocks N, we have from (13), 



C* = 6 [B log log 



K 
B 



(14) 



The following theorem gives similar results for other fading 
models. 

Theorem 5. If \vf k \ belongs to either Nakagami-m, Weibull, 
or LogNormal family of distributions, then, for regular 
extended networks, the scaling laws for the upper bounds 



For Nakagami-(m, w): C* = 0(BN log log -§) 

For Weibull(\, t): C* = 0(BN log log^ -§) 
ForLogNormal(a,uj): C* = O^BN^J\og§ 

and the scaling laws for the lower bounds are: 

For Nakagami-(m, w): C* = fl (BNf^ N (r, N) log log 



For Weibull(X, t): C* = O (BNf^ N (r, N) lo g log t f ] 
For LogNormal(a, u): C* = O ^BNf™ (r, N) y/log f 



Proof: Proof sketch is given in Appendix C. For detailed 
proof, see [16]. ■ 



Corollary 4. // » X> » i, and » 

1 for Nakagami, Weibull, and LogNormal distributions, re- 
spectively, then tighter upper bounds on the scaling laws can 
be found as follows: 

For Nakagami-(m, w): C* = O ^BN log 



For Weibull(\, t): C* = O (^BN log 



log t 



N 



For LogNormal(a, uS): C* = 0[ BN log cxp i^ s S_ l 



The proof sketch is given at the end of Appendix C. For 
detailed proof, see [16]. 

V. Design Principles 

In this section, we interpret the obtained results and outline 
some fundamental design principles for the service providers 
and network designers in order to achieve provisioning of 
various QoS guarantees for the end users and at the same time 
maximize the revenue for the service providers. For simplicity, 
we only consider Rayleigh fading channels, though, we note 
that similar results can be obtained for different fading 
models. We focus on four scenarios in this paper. In the 
sequel, we call our system scalable under a certain condition, 
if the condition is not violated as the number of users K — > oo. 

Principle 1. In dense femtocell deployments, with the 
condition that the per-user throughput remains above a 
certain lower bound, for the system to be scalable, BN must 
scale as ft( log g gJ J . 

We use the dense-network abstraction for a dense femtocell 
deployment [17] where the service operator wants to maintain 
a minimum throughput per user. In such cases, a necessary 
condition that the service provider must satisfy is: 



BN[ log ( 1 + Pcon/3 2 r-o za log 
K 



Kr 



> S, 



for some s > 0, where s = O(l). The above equation implies 

BN log log K 



= 0(1) 



(15) 



Therefore, the total number of independent resources BN, 
i.e., the product of number of base stations and the number 
of resource blocks (i.e., the bandwidth), must scale no slower 
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than i og ^ g K ■ Otherwise, then the system is not scalable 
and a minimum per-user throughput requirement cannot be 
maintained. 



Principle 2. In a large extended multi-cellular network, 
if the users are charged based on the number of bits they 
download and there is a unit cost for each base station 
incurred by the service provider, then there is a finite range 
of values for the user-density ^ in order to maximize return- 
on-investment of the service provider while maintaining a 
minimum per-user throughput. 

Consider a regular extended network with fixed number of 
resource blocks N. In this case, we have C* = Q(B log log ^) 
from (14). We assume a revenue model for the service provider 
wherein the service provider charges per bit provided to the 
users. Thus, the overall return on investment of the service 
provider is proportional to the achievable sum-rate per base- 
station. In large scale systems (large K), the associated opti- 
mization problem is: 



max c log 

K,B 



1 + Pcon/3 



-2m 



log 



s.t. 



cB log(l + P c 



2 

m /3 2 ro Za lo 



BR 2 ) 

BR 2 , 



K 



> s, 



(16) 



for some s > 0, where c is a constant bounded as described 
in (12)-(13). For simplicity of the analysis, let /3 = rrj = R = 
Peon = 1 and - = 0.1 (in respective SI units). Defining p = 
^, the above problem becomes a convex optimization problem 
in the variable p. Then, finding the optimal p via Lagrange 
multiplier method, the following Karush-Kuhn-Tucker (KKT) 
condition must be satisfied: 

(A + 1)10 



P = 



(17) 



(l + logp)A' 

where A > is the Lagrange multiplier. Note that, the La- 
grange multiplier represents the cost associated with violating 
the per-user throughput constraint. The plots of LHS and 
RHS of (17) along with the constraint curve are plotted for 
A = 0.1, 1, oo in Fig. 3. Here, the constraint curve (see the con- 
straint in (16)) is given by f log (1+ log p). Note that according 
to (16), the constraint is satisfied if p < S log(l + logp), i.e., 
when the constraint curve lies above the LHS curve. We notice 
from Fig. 3 that this occurs when p 6 [1.1, 12.7]. Since the 
optimal solution satisfies the constraint in (16), the optimal p 
lies in the set [1.1, 12.7]. Figure 3 also shows that as the cost 
of violating the constraint, i.e., A, increases, the optimal user- 
density p for a given A satisfying (17), i.e., the value of p at the 
intersection point of LHS and RHS curves in Fig. 3, decreases. 
Since p = 4.1 corresponds to A = oo, the optimal p is greater 
than or equal to 4.1. Summarizing the above observations, the 
optimal user density (number of users per base-station) lies in 
the closed set [4.1, 12.7], a finite range of values as stated in 
Principle 2. 

To investigate the variation of optimal p for a given A, 
denoted by p*(A), we plot p* (A)-versus-A in Fig. 4. As shown 
earlier, the optimal user density lies in a finite range (here, 



LHS 

RHS for X. 0.1 

-O- RHS for X - 1 
RHS for A = ~ 
-!>- Constrain! curve 




10 10 

User density p = K/B 



Fig. 3: LHS and RHS of (17) as a function of p. 




10 10 

Lagrange multiplier A 



Fig. 4: Optimal user-density, i.e, p*(X), as a function of A. 



between 4.1 and 12.7 users/BS for A > 0.29). Furthermore, 
the optimal user-density p* (A) is a strictly-decreasing convex 
function of the cost associated with violating the per-user 
throughput constraint, i.e., A. 

Principle 3. In a large extended multi-cellular network, 
if the users are charged a fixed amount regardless of the 
number of bits they download and there is a unit cost for 
each base station incurred by the service provider, then there 
is a finite range of values for the user-density in order to 
maximize return- on- investment of the service provider while 
maintaining a minimum per-user throughput. 

Consider a regular extended network with fixed number of 
resource blocks N, similar to that assumed in Principle 2. 
Here, we assume a revenue model for the service provider 
wherein the service provider charges each user a fixed amount 
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regardless of the number of bits the user downloads. Then, the 
return on investment of the service provider is proportional to 
the user-density p = ^. In large scale systems (large K), the 
associated optimization problem is: 

K 

max s — 

K,B B 



s.t. 



cSlog(l + P C on/3 2 r - log 



BR' 2 



K 



> 8 



(18) 



for some constants c, s, s > 0. Here, s depends on the amount 
users are charged by the service provider, and c can be 
bounded according to (12)-(13). For simplicity of analysis, 
let j3 = To = R = Peon = 1 (in respective SI units). 
Similar to Principle 2, the above problem becomes a convex 
optimization problem in the variable p = ^. Therefore, the 
optimal solution, denoted by p* , satisfies the constraint in (18) 
with equality. In particular, we must satisfy 



£ < log(l 



logp) 



(19) 



for all feasible values of p. The plot of LHS and RHS of (19) 
as a function of p (for p > 1) is plotted in Fig. 5. Examining 




A B 



10 10 

User density p = K/B 



Fig. 5: LHS and RHS of (19) as a function of p. 

(19) and Fig. 5, we note that the per-user throughput constraint 
is satisfied only if - G [0,0.26]. Moreover, for a given value 
of -, p lies in a closed set (for which the RHS curve remains 
above the LHS curve). The maximum value of p in this closed 
set, i.e., the value of p at point B in Fig. 5, is the one that 
maximizes the objective in (18), i.e., sK/B. Hence, it is the 
optimal p for the given value of s/c. Let us denote it by 
p*{s/c). Note that p*(s/c) > 2.14 (since point B lies to the 
right of point A in Fig. 5). 

If s/c is known exactly, then the optimal user-density 
p* = p*(s/c). If not, we can write from (12)-(13) that 
cib < c < c u b, for some positive constants cib, c u b- Then, 
p* £ [p*(s/c\ b ), p*(s/c u b)}. Moreover, since p*(s/c) > 2.14 



for all s/c e [0, 0.26], we have p* (s/c ub ) > p* (s/c\ b ) > 2.14. 

Principle 4. In regular extended networks, if the users 
are charged based on the number of bits they download and 
there is a unit cost for each base station and a cost cn for 
unit resource block incurred by the service provider, with the 
condition that return-on-investment remains above a certain 
lower bound, then for fixed B, the system is scalable only if 
N = 0(logK), and for fixed N, the system is scalable only 
ifB = 0(K). 

Consider the case of a regular extended networks with large 
K. Using (8) in conjunction with (40) (using the upper bound 
obtained via Jensen's inequality in (8)), we have 



C* < log 1 



Pc, 



N 



Ik + 



P a 



N 



| +0(1) jBN 



BN log (^l K 



for large 



Pconl 



con'K 



(20) 



where Ik = fi 2 rQ 2a log BI ^° . For simplicity of analysis, let 
^con = /3 = rrj = -R = l(in their respective SI units). If 
the service provider wants to maintain a minimum level of 
return-on-investment, then 



BN 



i ( 1 i KN \ - 



B + c N N 

for some s > 0. The above equation implies 

N = O(logA'), for fixed B, and 
B = O(K), for fixed N. 

VI. Conclusion 



(21) 



(22) 



In this paper, we developed bounds for the achievable 
downlink rate in multi-cellular OFDMA based networks and 
specified the associated scaling laws with respect to number 
of users K, number of base-stations B, and number of 
resource-blocks N (bandwidth). Our general bounds hold for a 
general spatial distribution of base-stations, a truncated path- 
loss model, and a general channel-fading model. We evalu- 
ated the bounds for dense and extended networks in which 
nodes were distributed uniformly for Rayleigh, Nakagami- 
m, Weibull, and LogNormal fading models. We showed that 
for dense networks, under Rayleigh fading channels, the 
achievable rate is lower bounded by Sl(min(i?, N) log log K), 
and upper bounded by 0(BN log log K). The corresponding 
result for regular extended networks showed that the capacity 
is lower bounded by Q (Slog log and upper bounded by 
0(BN ^loglog We derived similar results for Nakagami- 
to, Weibull, and LogNormal family of fading models. 

We then applied the obtained results to develop four design 
principles for the service providers and regulators to achieve 
QoS provisioning along with system scalability. According to 
the first principle, in dense-femtocell deployments, if a min- 
imum per-user throughput requirement must be maintained, 
then the system is scalable only if BN scales as Q ( log ^ K ) . 
In the second and the third principles, we considered two user- 
charging methods: per-bit and fixed, and showed that the user 
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density must be kept within a finite range of values in order 
to maximize the return on investment, while keeping the per- 
user rate above a certain value. Finally, in the fourth principle, 
we also considered the cost of the bandwidth to the service 
provider along with the cost of the base stations and showed 
that for fixed B, the system is scalable only if N = 0(log K), 
and for fixed N, the system is scalable only if B = 0(K). 

Appendix A 
Proof of Theorem 1 

By ignoring the interference, we have 

B N 



i—1 n—1 



B N 



have Efc ^k( x ' 2/, t^) = 1 V i, n. Now, (29) can be re-written 



C x , y ,AU,P) < EE l0 S( 1 + ^^(n)J 

This implies 

maxE 

{ log (1 + 7&)} 

i—1 n—1 
B N 

^EE E { m f lo s( 1 + Pcon7,"fe) } (24) 

i—1 n—1 
B N 

<^^E{log(l + P con max 7 £ fe ) } , (25) 

i—1 n—1 

where (24) follows because, for any function /(■,•), 
maxfe E{/(fc, •)} < Ejmaxfe f(k, •)}, and (25) follows be- 
cause log(-) is a non-decreasing function. One can also con- 
struct an alternate upper bound by applying Jensen's inequality 
to the RHS of (23) as follows: 

B 



1=1 

B 



N 

n 

Peon 



< iV^log l + ^max 7 ^), (27) 



N n,k 



since En -P™ - Pon- Therefore, 

B , 

c* < ivE E 1 lo %{ 1 

1=1 



Pc< 



(28) 



Combining (25) and (28), we obtain the upper bound in 
Theorem 1. 

For lower bound, let P C on/N power be allocated to each 
resource-block by every BS. Then, 



Cx,y,v(U . P) 



B N 



>EE lo « 1 



-fcon 7; 



ki(r 



N 



Peon Ej^i Ylk, 



(29) 



i=l n=l "' ' --""^3^1 'J,fti(n)- 

where ki(n) is any other user allocated on subchannel n by 
BS i. Note that, due to sub-optimal power allocation, all user- 
allocation strategies {ki(n), Vi, n} achieve a utility that is 
lower that C x , y ,„(U, P). To handle (29) easily, we introduce 
an indicator variable I™ k (x, y, v) which equals 1 if k = ki(n), 
otherwise takes the value 0. Since, each BS i can schedule at- 
most one user on any resource block n in a given time-slot, we 



x,y,v(U , P) 

B N K 



(30) 



^ EEE^^y^) 1 ^ 1 



P <v" 

-fcon 7j k 



N- 



Taking expectation w.r.t. (x,y,is), we get 

P=on7" fc 



C* > E I I^ k (x,y, z)log(l 

i,n.k \ 



N + Peon Ei#i 7 



>VBjj ?1 (.,,,,)^i±^l-V (31) 

Here, the last equation holds because for any non-decreasing 
concave function V(-) (for example, V(x) = log(l + x)) and 
for all di, d,2 > 0, we have 

'd x - 



Now, 



JV 



^y&)> v ™r m + v(o). 

— t- < 1. Therefore, 



c* > E E ' 



' JT^s, log (1 + Peon 7 fr 
^ + PoonE jW 7^ fe 



(32) 



(33) 



To obtain the best lower bound, we now select the user ki (n) 
to be the one for which 7 ™ fe attains the highest value for every 
combination (i, n), i.e., 

IU*^) = \\ ** = (34) 
I otherwise. 

Using (34) in (33), we get the lower bound in Theorem 1. 

Appendix B 
Proof of Theorem 2 and Theorem 4 

In this proof, Lemma 1, Lemma 2, Lemma 3, and Lemma 4 
have been stated without proofs. For detailed proofs, see the 
technical report [16]. 

To prove the final result, we will use three lemmas. 
Lemma 1 gives a new lower bound from (8) using an ap- 
plication of one-sided variant of Chebyshev's inequality (also 
called Cantelli's inequality). 

Lemma 1. In a dense-network and Rayleigh fading channel, 
i.e., fT*. ~ CAf(0, 1), the lower bound on the achievable sum- 
rate of the system is 

C* (35) 
> f™(r, B, N) £ E x , y , u { log (1 + P con max 7 ^)}, 

z,n 

where r > is a fixed number, and f® N (r, B, N) = 

^— — — - — r „„ _•,„, r-^r, and u = a = 1 for Rayleigh 

fading channels. 
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Lemma 1 in conjunction with Theorem 1 shows that both, 
the upper and lower bounds, on the achievable sum-rate, i.e., 
C* , are functions of SNR scaling. To find the SNR-scaling 
law, we state the following lemmas. 

Lemma 2. Under Rayleigh fading, i.e., v™ k ~ £/V"(0, 1), the 
CDF of 7,™ fe is 



n2 



af3 2 p 2 



exp(-j/g)ds(g), 



(36) 



where d = y/aj + b 2 , and s(g) equals 




Using Lemma 2, we can now determine the scaling of 
maxfe7™ fe for a given (i,n) under Rayleigh fading. 

Lemma 3. Let 7T 1 , be a random variable with a cdf defined 



in (36). Then, the growth function 



i--F T » (7) 



converges to a 



constant f3 2 r 2a as 7 — > 00, and 7 ™ fc belongs to a domain 
of attraction [15]. Furthermore, the cdf of (maxfc j" k — Ik) 
converges in distribution to a limiting random variable with a 
Gumbel type cdf, that is given by 



cxp(-e xr o a /P 2 ^ x e (-00,00), 



(37) 



where Ik is such that Fyn (Ik) = 1 — 1/A. In particular, 
l K =p 2 r^\og^. 

Note that the SNR scaling is independent of (di,bi) and 
depends on ro,p. Now, we use the above three lemmas to 
prove the final result in Theorem 2 (and Theorem 4). Since 
the growth function converges to a constant (see Lemma 3), 
we apply [9, Theorem A. 2] giving us: 



Pt{Ik ~ log logiC < max7™ fc < l K + log log A"} 

k ' 



> l-O 



1 

logAry 



(38) 



where l K = /3 2 r^ 2a log Now, 



E { log (1 + P C on max7&)} 
< Pr{ max 7™,. <l K + log log A | 
x log(l + P C0 Jk + Peon log log A) 
+ Prj max7™ A . > l K + log log A" | 
x log(l + Pcon/S 2 ^ 2 " A) 



(39) 



< log(l + Peon I K + Pcon log log A") 

+ o( I ^)log(l + P CO n/3 2 r o : 2Q A) 

< log(l + P con l K + P con loglog A) + 0(1), (40) 

where, in (39), we have used the fact that the sum-rate is 
bounded above by log(l + P con (3 2 rQ 2a K). This is because 

E { log (1 + P con max 7^)} 

< E { log (l + Peon 

k 

< log (l + P con 5>{ 7 £j) 

k 

< log(l + P con /3 2 r - 2Q ^E{|< fc | 2 }) 

(41) 



< log (l + Peon/3 2 ro 2Q A 
Further, from (38), we have 
E { log (1 + P^ max 7^)} 

> log(l + P con l K - Peon log log A) (l-O 



(42) 



1 



log A 



Combining (40) and (42), we have for large A, 
log(l + PcoJk) + 0(1) < E { log (1 + P con max 7^,) } 
< log(l + P C0n Z K ) + O(l), (43) 

Therefore, 

E{ log (1 + Pcon maxfcT^)} 
log log A 

By definition of G(-) in Section I, we have, 

E { log (1 + P con max 7 £ fc ) } = 6 ( log log A) . (45) 

Using the above equation in conjunction with (8) and (35), we 
have 



lim 

K — >oo 



(44) 



C* = 0{BN log log A), and 



C* = n(BNf° N (r : B : N)\oglogK), 



(46) 



where f^ N (r, B, N) is defined in Lemma 1. 

Now, to prove Corollary 1 and Corollary 3, we use the upper 
bound in (8) obtained via Jensen's inequality. In particular, we 
have 



C* < N 



E e (^('+t?*))| (47) 
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lim 



This implies 



E {log(l + Sf max r ^ 7 ^)} 

1 at 1 ■ 



, log KN 



C* = 0(BN\og 



log KN 
~N 



(48) 



(49) 



Note that the above result is only true if P con log ^s> 1 (or 
if, log £jtf- 3> 1 since P CO n is fixed). 

To prove Theorem 4 for regular extended networks, we use 
the following lemma instead of Lemma 1. 

Lemma 4. In a regular extended-network and Rayleigh fading 
channel, i.e., v™ k ~ £/V"(0, 1), the lower bound on achievable 
sum-rate of the system is 

C* > J lo N (r, N) J2 E { log (1 + Pon max 7 £ fc ) } , (50) 



where r > is a fixed number, f™(r, N) 



p o2 2-2, 

c = (4 



N+(/i+rcr)c 



, and 



R 2 V" ' v / 3(2a-2) / 

The rest of the steps in the proof remain same, with the 
only change of substituting p 2 by R 2 B in Lemma 3 and in 
the subsequent steps (38)-(49). 

Appendix C 
Proof Sketch of Theorem 3 and Theorem 5 

For complete proof, see the technical report [16]. The proof 
comprises of three parts. First, we take the limit of the growth 

l-*>, (7) 



function /i( 7 ) 



as 7 — > 00 and show that the 



distribution of 7™^ lies in the domain of maximal attraction 
[15]. In particular, we determine if the limiting distribution of 
7™ fc is a Frechet, Weibull, or Gumbel type, and find Ik in each 
case such that 1 — F^? {Ik) = 7^- Second, using the limiting 
pdf (Frechet, Weibull, or Gumbel), we show that max^ 7™ fc = 
0(Ik) with high probability. This can be proved by following 
the steps of [9, Theorem 1] for Gumbel type distributions. 
Finally, we use the same steps as used (39)-(46) to find the 
scaling laws of C* . We briefly give the main calculations, for 
Nakagami-m, Weibull, and LogNormal distributions here. 



1) For Nakagami-(m, w) fading, h{j) 



^/3V as 

7 —t 00. Hence, /i'( 7 ) — > and the limiting dis- 
tribution of max/j 7™ fe is of Gumbel type. Solving 
for Ik, it can be shows that for large K, Ik = 

2) For Weibull-(A, t) fading with scale A > and shape t > 

0, h(n) = 2lX '£ t r ,C? t/2 ■ Therefore, lim^ h'(j) - 
0. Consequently, the limiting distribution of maxj 7"^ 
is of Gumbel type. In this case, 

/ A - = A 2 /3 2 r - 2 "log?^. 

3) For LogNormal fading, i.e, log Af(a, w), we have ^(7) w 

for large 7. Therefore, lim 7 _ i . ocl h! (7) = 0, and the 



limiting distribution of is of Gumbel type. In this 
case, we finally get 



K 



Slj log 



L+e(io g io g K) 



Note that p is fixed for dense networks and p m r\[~B 
for regular extended networks. The proofs of Corollary 2 
and Corollary 4 follow from the proofs of Corollary 1 
and Corollary 3 in (47) - (49) in Appendix B (by using 
the upper bound on the achievable sum-rate obtained via 
Jensen's inequality in Theorem 1). 
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Detailed proofs 



I. Proof of Theorem 2 and Corollary 1 for Dense Networks 

The proof outline is as follows. We first prove three lemmas. The first lemma, i.e, Lemma 1.1, uses one-sided 
variant of Chebyshev's inequality (also called Cantelli's inequality) and Theorem 1 to show that 



C > f° N (r,B,N)J2v{ log (1 + Peon max 7 £ fc )}, 



where C* is expected achievable sum-rate of the system. The second lemma, i.e, Lemma 1.2, finds the cumulative 
distribution function (CDF) of channel-SNR, denoted by F 7 » (■), under Rayleigh-distributed and a trun- 

cated path-loss model. The third lemma, i.e, Lemma 1.3, uses Lemma 1.2 and extreme-value theory to show that 
(maxj; Yik ~ Ik) converges in distribution to a limiting random variable with a Gumbel type cdf, that is given by 

cxp^e""""^ 2 ), x G (-00,00), (1.1) 

where F^n {Ik) = 1 — Thereafter, we use Theorem 1, Lemma 1.1, Lemma 1.3, and [SH05, Theorem A.2] to 
obtain the final result. 

Now, we give details of the full proof. 

Lemma 1.1. In a dense -network, the expected achievable sum-rate is lower bounded as: 

C* > /r(r, J B,iV)^E{log(l + P con max 7l l ; fe )} ! (1.2) 

i,n 

where r > is a fixed number, fP N (r, B, N) = -, jt— — - — r _,„ _ irt , — — , 11 and a are the mean and standard- 



deviation of |i/ n 



(l+r 2 )(iV+i 3 C o„ l 9 2 ro^ a (|i+r<7)B)' 



Proof: We know that 

E^ = ^ 2 E *s2>j:*l a < W a E K^i 2 - cl3) 



Therefore, the lower bound in Theorem 1 reduces to the following equation. 

max fc log ( 1 + Peon 7"fc ) 1 



^>E E { 

i,n,k \ 



N + P C on/? 2 r - 2a 



n 12 



(1.4) 



n 12 



Now, we apply one-sided variant of Chebyshev's inequality (also called Cantelli's inequality) to the term \ v j k 

k 



in the denominator. By assumption, |^™ fc | 2 are i.i.d. across i,k,n with mean // and variance a. Hence, applying 
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Cantelli's inequality, We have 

=► Pr(Y,Kk\ 2 >(» + r<j)B) < — Lj (1.5) 
=> Pr(j2Kk\ 2 <(^ + ra)B) > ^ (1.6) 



where r > is a fixed number. 

Now, we break the expectation in (1.4) into two parts — one with Yljjti Wj'kl 2 > (m + r<r)B and other with 
Tlj^i \ vl jk\ 2 — + r<r)B. We then ignore the first part to obtain another lower bound. Therefore, we now have 

B N 

N + (» + ra)BP con p2ro 

B N 



z— 1 n— 1 



> 



1+r 2 



o- VVe( max log ( 1 + P con Y,\. )\ (1.7) 

= /,™(r> B,N)J2^{ log (1 + ^con max^)}, (1.8) 

i,n 

where (1.7) follows because Y^jjti Wj'kl 2 ' s independent of vf h (and hence, independent of 7™ fe ). Note that for 

Rayleigh fading channels, fj, = a = 1. ■ 
Lemma 1.1 and earlier proved Theorem 1 show that the lower and upper bounds on C* are functions of max/. 7™ fc . 
To compute maxfc 7™ ft for large K, we prove two lemmas, Lemma 2 and Lemma 3. 

Lemma 1.2. Under Rayleigh fading, i.e., v™ k ~ CA/ r (0, 1), //ze CDF of is given by 



i e 2 .. 

-x ( g , f <j>-w* , 

(p+ti) 2 " 



^,(7) = l -f2^^ -^ifll * ' / - exp(- 7 / 5 )d s ( 5 ), (1.9) 



where d = \/ af + b 2 , and 



<9) = — 2 



1 



c « s ^—77^ +p 2 cos- 1 



(JO ~" V 2d(-^)- 1/2Q y ' " ~" V 2dp 



Proof: We assume that the users are distributed uniformly in a circular area of radius p and there are B 
base-stations in that area as shown in Fig. 1.1. 

The probability density function of the user-coordinates {xk,Vk) can be written as 

!-^ x 2 + y 2 <p 2 
(LIO) 
otherwise. 
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Note that around any base-station, the users are distributed at-least within a distance R (R > tq). Hence, p — d 
P ~ V a f + > R > r o f° r a U i- Now, 



Ilk 



i{r 0)V /(a; fe -a i ) 2 + (y fe -& i ) 2 }) Vl^l 



(1.11) 



(1.12) 



We now compute the probability density function of Gi t k (= Z^ 2 ^ l"). 



-l/2a\ 



if 9 > P 2 ro 2a 



(1-13) 



otherwise. 



Pr( v /(x fe -a l ) 2 + (y fe -6 J ) 2 < (£) V2Q ) 

Now, Pr^^f^'fc — a i) 2 + ~ < (^) 1//2a J is basically the probability that the distance between the user 
k and BS i is less than (jfi) 1 ^ 2a . Since, the users are uniformly distributed, this probability is precisely equal to 
^-5 times the intersection area of the overall area (of radius p around O) and a circle around BS i with a radius 
°f (^f?) 1 ^ 2q . This is shown as the shaded region in Fig. 1.2. 
Therefore, we have: 



Pr(G,, fe > g) = { 



1 


if(f) 


-1/2q 


€ (p + ef, oo) 


s(g) 




-1/2q 


e (p — d, p + d] 


(#) 




-1/2q 


e (r ,p- cf] 







-1/2q 


G [0,r ], 



(1.14) 
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Fig. 1.2: System Layout. The BS i is located at a distance of c? from the center with the coordinates (a,i,bi), and 
the user is stationed at {xk,Uk)- 



where 



s(g) = — j 

7Tjr 



V/3 2 



-l/c 



COS 



' 3 \- 1 / c 



-l/2a 



1 — 1 
■p COS 



d 2 +p 2 



2dp 



-l/2a 



-l/2a 

4t) -d)(d+(-^ 



/3 2 



-l/2o 



The CDF of G^k can now be written as 

if .9 e [0,/3 2 (p + d)- 2a ) 
l-s(g) ifge [(3 2 (p + d)~ 2 <*,(3 2 (p~d)- 2a ) 
i-(^r 1/a ^ if ge[p 2 ( P -dr 2 «,p 2 r^) 

1 if 5 e [/3 2 r - 2Q ,oo) ; 
A plot of the above CDF is shown in Fig. 1.3. 

The probability density function of Gi.k can be written as follows: 

if g€[0,l3 2 (p + d)- 2a ) 

ifge[p 2 (p + d)- 2 <*,P 2 (p-d)- 2 >*) 
^(^r 1 - 17 " ifg€[^( P -d)-^,^) 

if g = (3 2 r^ 2a 
ifg>P 2 r- 2a 1 



p 2 



(1.15) 



-l/2a\ 



(1.16) 



(1.17) 



where ds j^ < 0. The pdf of Gi.k has a discontinuity of the first-kind at /3 2 r 2a (where it takes an impulse value), 
and is continuous in [f3 2 (p + d)~ 2a , /3 2 r^~ 2a ). At all other points, it takes the value 0. 
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Pr(G hk < g) 



1 




Fig. 1.3: Cumulative distribution function of G^fc. 



Using (1.17), the cumulative distribution function of 7™ fe , i.e., F 7 n (7) (when 7 > 0) can be written as 

*V&(7) 

= /p(l< fe | 2 <T)f GiJt {g)dg (1.18) 
(l- e -7/s)/ G . fe ( 5 ) dff (1.19) 

e- j/9 f G U9)dg (1.20) 



1 - %e ^ - \ e -"»-±s{±] f'-^dg + / e^'d»(g). (1.21) 



l--F 7 « (7) 

Lemma 1.3. Let 7™ ft fee a random variable with a cdf defined in Lemma 1.2. Then, the growth function j n ''^ 
converges to a constant /3 2 rQ 2a as 7 — > 00, and 7™ fc belongs to a domain of attraction [DN03]. Furthermore, the 
cdf of (maxfc 7™ fe — Ik) converges in distribution to a limiting random variable with a Gumbel type cdf, that is 
given by 

cxp(-e- xr ° a/l32 ), x G (-00,00), (1.22) 
where Ik is such that [Ik) = 1 — 1/K. In particular, Ik — 1 r^ ) ' la log -^-J 2 -. 
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Proof: We have from Lemma 1.2 

F y?M = 1 



r 2 7 



1 , g s-i-i/ 



^(p-d)- 2 "" a/3 2 p 2K /3 2 



P 2 (p-d)- 



P 2 (P+dY 



e-i/°s'(g)dg 



1 



r 2 7 



fl 2 — 2o 



/3 2 (p-d) 



e-^^iir^dg + s^e-^ 



/3 2 (p+d)- 2 ' 
r 2 7 

9 



— 2 ^ 



e s ■ 



1 ,Sr'- 



/3 2 (p-d)- 
/3 2 (p+d)- 



(1.23) 



p2 {p _ d) -2 a ~ af3 2 p 2X f3 2 

<£j_s{g) 

p 2 { P +d)- 2 <* 9 



dg, 



{p-df 



where < ^ < s(g) < I (see Fig. 1.3). Now, we claim that 



lim (l — F~n h))e^ 2r o- a = 



2 

2 ' 



(1.24) 



(1.25) 



It is clear that the first two terms in (1.24) contribute everything to the limit in (1.25). We will consider the rest of 
the terms now and show that they contribute zero towards the limit in RHS of (1.25). First, considering the 4 th , 5 th , 
and 6 th terms, we have 



lim e* 9 r o x 

7— >-oo 



7 



< lim i 

7— ¥ OO 



g /3 2 (p-d)- 



{P-df 



g /3 2 (p + d)- 



P 2 (p+d)- 2a 9 2 

7 



dg 



e °s(g) 
7fi dg 



< l im (P~f e -^( P -^-^ +e -- 
= 0. 



_ 7 ((p-rf) 2 "-r 2 °) 
■((p+d) -r ) e 



/? 4 (p + rf) 



-4a 



(1.26) 

(1.27) 

(1.28) 
(1.29) 



Now, we consider the third term in (1.24). We will show that 



lim e ,3 ' lr o " x 

7— >oo 



1 f 9 



d 9 = 0. 



Taking the exponential inside the integral, we have 



T( 7 ) 

Substituting 7/5 by x, we get 

T(7) = 



2 r- 2 ° 



/3 2 (p-d)- 



T(7) 



= -x+7r 2 °//3 2 



7 



af3 2 p 2 \xP 2 



-l-i 



— tt \ dx. 



(1.30) 



(1.31) 



(1.32) 
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Again substituting x ~ jr^ //? by y, we have 



ap 2 \p 2 J J V p 2 

( (p-d) 2a -r$<>)-, 

-i 1 2r> i _i_ 1 — — 

(j_Y°(inL^~ - 



P 2 „ . „/ ..(P-<J) 2 °-'-n 



-r - 2Q+2 1-e" 7 ^ ), (1.35) 

where, in (1.34), an upper bound is taken by putting y = in the term (a) + inside the integral. Since 

T(7) is positive, (1.35) shows that lim 7 _ i . 00 T("f) = 0. Hence, the claim is true. 

Now, after computing the derivative of -F 7 « (7) w.r.t. 7 to obtain the probability density function f 1 ^ k (7), we 
have 

Urn f^)eri a lf = r l^-. (1.36) 

We do not prove the above equation here as (1.36) is straightforward to verify (similar to the steps taken to prove 
(1.25)). From (1.25) and (1.36), we obtain that the growth function converges to a constant, i.e., 

1 - F„n (7) 

lim , 7 ': fc =P 2 r^ a . (1.37) 

The above equation implies that 7™ fc belongs to a domain of maximal attraction [DN03, pp. 296]. In particular, the 
cdf of (maxfe7™ A . — Ik) converges in distribution to a limiting random variable with an extreme-value cdf, that is 
given by 

cxpt-e-"^ 2 ), x G (-00,00). (1.38) 
Here, Ik is such that (l K ) = 1 — 1/K. Solving for Ik, we have 

> + /t «-* 5 M£)- , - i *+/f V *<-'<»>>* (,J9) 



K p 2 



(p-<i) 2 
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Substituting Ik 1 9 by x in the first integral in RHS of (1.39) and computing an upper bound, we get 



1 

K 



< 'JL e 



a/3 2 p 2 



[k \xB 2 ) \ x 2 J 



dx - e P 2 ip- d r 



(p + d) 2 °> 



s'{g)dg (1.40) 



exp 



/ Ik 

V B 2 r^ 



1 /I 



B 2 ro 2a Jp 2 ' ap 2 \p 



K 



-is 



e x x 1+ idx + e e 2 (p- d )- 



' 



(p-dH 

(P+3T 2 



< -4 



L e 2 

i2 



ap 2 V /3 2 / 



2a\ 



V P 2 



e x dx + e ' 32 (f- d >" 



< e 



r5« r o 



„2-2a 



< e 



p 2 ap 2 V /3 



(a 



e 31 da: + e 2 (p-<*)- 



m 



e 



r- 2 



ap 2 V /3 



1 + 



a/it 
1 

Ik~ 



y e /3 2 (p-d) 



e 

,,2 



1.41) 

(1.42) 

(1.43) 
(1.44) 
(1.45) 



In (1.40), we substitute Ik 1 9 by x in the first integral of (1.39), and compute an upper bound by taking the exponential 
term out of the second integral of (1.39). In (1.41), we note that ^ p ~f* < s(g) < 1. From the above analysis, we 



now have 



Kr 2 

2 -2a i A "0 



P 



Now, to compute a lower bound on from (1.39), we note that fact that ds ,^ < 0. Therefore, 



r 2 <K 

R p z 



Kr 2 



Ik > B 2 r^ 2a log 



From (1.46) and (1.48), we have 8 2 r„ 2a log ^<l K < P 2 ^ 2 " log ^ + O(i^). Therefore, 



(1.46) 

(1.47) 
(1.48) 

(1.49) 



for large K. ■ 
Interestingly, for a given BS i, the scaling of maxj 7™ fe (given by Ik in large AT regime) is independent of the 
coordinates (dj, 6,) and is a function of ro,p. Now, since the growth function converges to a constant (Lemma 1.3), 
we apply [SH05, Theorem A. 2] giving us: 



Pr[l K - log log A' < max 7^ < l K + log log k] > 1 - o(j^) , 



(1.50) 
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where Ik = (3 2 r Q 2a log Therefore, 



E { log (1 + Peon max 7^)} < Pr(max7™ fc < l K + loglogif) log(l + PWk + P con loglogiT) 



+ Pr ( max 7 £ fc > l K + log log Kj log(l + P C0 ^ 2 r^ 2a K) (1.51) 

< log(l + P con l K + Peon loglogX) + log(l + P CO n/3 2 r - 2Q A') x o(^^r) 

= log(l + Pcon^A") + 0(1). (1.52) 

where, in (1.51), we have used the fact that the sum-rate is bounded above by log(l + Pcon/3 2 ?"n~ 2c * K). This is 
because 

E { log (1 + P con max 7 £ fe ) } < E { log (l + P con £ 7^) } 

fc 

< log (l + P con J2 E ^i,k}) (1-53) 

k 

< log^ + Pcon^^Ell^l 2 }) 

k 

< log(l + P con p 2 ro 2a K^), (1.54) 
where (1.53) follows from Jensen's inequality. Further, from (1.50), we have 

E { log (1 + P con max 7 £ fc ) } > log(l + P con l K - P con log log K)(l- O (—^j ) • (1-55) 
Combining (1.52) and (1.55), we get, for large K, 

BN log(l + P con ^- Pcon log log K){l-0{^^j) < ^E{ log (1 + P^ max 7^)} (1.56) 

< (log(l+P C0n l K )+O(l))BN. 

Therefore, from Lemma 1.1 and Theorem 1, we get 

(\og(l + P con l K ) + 0(l))BNfg N (r,B,N) <C* < (log(l + P con l K ) + 0(1)) BN (1.57) 
This results in: 

C* = 0(BN log log K), and 

C* = Q(BNf^(r,B, N)loglogK). (1.58) 

Now, to prove Corollary 1, we use the upper bound in Theorem 1 obtained via Jensen's inequality. In particular, 
we have 

C* < tf^E{bg(l + ^max 7 &)}} (1.59) 
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L9 

where (1.60) follows from (1.56), and Ikn = /3 2 r^~ 2a log K ^ 2 r " determines the SNR scaling of the maximum over 
KN i.i.d. random variables. This implies 

'* = C)( R/VW — 

N 

Note that the above result is only true if ^j^lxN ^S> L (or, log ^ff- 1 to make the approximation log(l + x) ps 
log a; valid for large x). 

II. Proof of Theorem 3 and Corollary 2 for dense networks 

We will first find the SNR scaling laws for each of the three families of distributions — Nakagami-m, Weibull, and 
LogNormal. This involves deriving the domain of attraction of channel-SNR 7™ fc for all three types of distributions. 
The domains of attraction are of three types - Frechet, Weibull, and Gumbel. Let the growth function be defined as 

l-F^n (7) 

h{"f) = j „''( fc 7 ) ■ The random variable, y? k , belongs to the Gumbel-type if lim 7 _ i . 00 ^'(7) = 0. It turns out that 
all three distributions considered in [AKS12], i.e., Nakagami-m, Weibull, and LogNormal, belong to this category. 
After showing this, we find the scaling, Ik, such that F 7 » (Ik) = 1 — 1/K. The intuition behind this choice of 
l K is that the cdf of max fc 7™ fc is Ffi (7). For 7 = l K , we have Fj^ (Ik) = (1 - ^/K) K -> e _1 . The fact that 
F^n (7) converges for a particular choice of 7 gives information about the asymptotic behavior of maxt 7™ fc . 

A. Nakagami-m 

In this case, \v£A is distributed according to Nakagami-(m. w) distribution. Hence, \v™ k \ 2 is distributed according 
to Gamma- (m, w/m) distribution. The cumulative distribution function of 7™ fe , i.e., F 7 r» (7) (when 7 > 0) is 

^ fe (7) = Jp(K k \ 2 <^)fG,. k (g)dg (Hi) 



r(m) 



fG i<k (9)d9 (H.2) 



/•Z^ 20 r(m, 21) 

1-/ f G ^ )dg (IL3) 

i/3 2 (p+d)- 2 ° 1 ( TO J 



where /c i ,.(5) is defined in (1.17). Now, for large 7, we can approximate (II. 3) as 



/3 2 r — 2a 

F 7r fc (7) * ^ f??) m ~ 1 e"-/^(5)^ 01.4) 



1 - r 2 / m 7 _ _J_ /to7\™~i ^ 1 ,_g^-i-A 

p 2 T(m) V w/3 2 r -2^ r(m ) Jp 2(p _ d) _ 2a \ wg J app? V 2 j 



o2 -2t> 



1 fP 2 (r-d)- 2a rm^m-i 



frrvyy 



'ds(g), (II.5) 
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where fa t k {g) is defined in (1.17). We claim that 



lim 7 1 - m e™ 32 -o- 2a 



e ™ 9 fG,, k (g) d 9 



(11.6) 



(11.7) 



p 2 T{m){wP 2 rQ 2a ) m - 1 ' 

Note that the first two terms in the RHS of (II. 5) contribute everything towards the limit in (II. 7). We will show 
that the rest of the terms contribute zero to the limit in RHS of (II.7). In particular, ignoring the constant T(m), 
the contribution of the two integral-terms (in (II. 5)) is 

P 2r a 2C " fm.rv\m-l m , 1 / n \ -1-4: /-/3 2 (p-d)~ 2 ° , mrj s m -l 



P 2 (p _ d) -2 a \wg / a/3 2 p 2 \p 2 



2 ro 2 ° 



f m \ 



m-l -IHJ.il. L. 

w ^ 9 H 2 T~ 



1 f a 



Jp 2 (p+d)- 2 " Kw 9 J 



Ti(7)+T 2 (7). 



Ti( 7 ) 



r 2 ( 7 ) 



ds(g) 



(II.8) 



Now, 



|Ti( 7 )| 



2 fl" 

m\ m - 1 /3« 



«; / ap* 



< 



(p-d) 2 ° 



g- m ~-e 



ix 1 1 l_ 



dg 



w / ap^ Jp-irl 



w / ap^ 

m x m-l 



(H.9) 
(11.10) 



J ^ max It) -br) 



/3- 2 r 2 ° 



0, as 7 — > oo 



2 m£1X 1 V £2 



,(p_ rf )2a, to+A-2 /r^ Q \'™+i-2] 1 - e-^l* 3 2 (p-d) 2 °-,3 2 r 2 = ) 



(£) 



(11.11) 

(11.12) 



where, in (11.10), we substituted i by x. Further, 



|T 2 ( 7 )I = 



6 



ds(g) 



(-^ ds(g] 



/3 2 ( P +d)- 2 ° \wg/ 



0, as 7 — > oo. 



(11.13) 

(11.14) 
(11.15) 
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Therefore, 71(7) ar, d 7a (7) have zero contribution to the RHS in (II. 7), and the our claim is true. Now, from (II.5), 
we have 

AfcW 01.16) 

m-l rP 2 r- 2 " fm yn _^ ( m _ 1)^-2 rP 2 ^ / m _ZZ . ... 

= pTT / ( — e lU9 fG t , k (g)dg — / — e <»<> f G (g)dg 

Using (II.6)-(II.7), it is easy to verify that 



lim /„« ( 7 ) 7 1 - m e™' 32 '-o- 2Q = ^ . (11.17) 

From (II.7) and (11.17), we obtain that the growth function converges to a constant. In particular, 

lim Y = , (11.18) 

7^00 / % „ fc ( 7 ) m 

Hence, 7™ fc belongs to the Gumbel-type [DN03] and maxfc 7™ fe — Ik converges in distribution to a limiting random 
variable with a Gumbel-type cdf, that is given by 

cxp(-e^ xr ° a/p2 ), x G (-00,00), (11.19) 



2 -2a 

- IrMT 

p 2 r(m)(w/3 2 ro 2 °)>- 

Now, since the growth function converges to a constant and Ik = (log If), we can use [SH05, Theorem 1] to 



where 1 - F 7 » fc (l K ) = -g. From (II.7), we have l K « ^ log p 2r( m v^ r -'°^-i for lar § e K - 



obtain: 



Pr I l K -log log X < max 7 ™ fc < Z K + log log if ) > 1 - O ( — !— 
I fe ' J VlogA 

This is the same as (1.50). Thus, following the same analysis as in (I.51)-(1.61), we get 



(11.20) 



C* = 0(BNloglog—£) and (11.21) 

C* = BNf° N (r,B,N)n(\oglog^y (11.22) 

Further, if log ^ > 1, then C* = O ( log ^^f j . 
B. Wej&w/Z 

In this case, |^™ fc | is distributed according to Weibull-(A, t) distribution. Hence, |^™ fe | 2 is distributed according 

to Weibull-(A 2 , t/2) distribution. We start with finding the cumulative distribution function of 7 ™ fc , i.e., F^n (7) 
(when 7 > 0) as 



(11.23) 



1-/ eTOT f Gik (g)dg (11.24) 

7^2( p+(i )-2 Q 



l--£e 



p 2 / f> 2 a(3 2 p 2 V/3 : 



(JO 1 C "dg+ f^ )a e (^) <k(0). (11.25) 
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This case is similar to the Rayleigh distribution scenario in (1.21). Therefore, it is easy to verify that 

C i ) t/2 2 
lim (1-F~n (j))e K,:,2 ''o' Aax2> = and (11.26) 

lim f^ k {l)l 1 - t,2 ^^^ ) = ~ f° W* ■ (IL27) 

l-F y n (7) 

Thus, the growth function h(j) = — f — '-f . can be approximated for large 7 as 

/i( 7 ) « KH °— '- — 7 1 -*/ 2 . (11.28) 

Since lim 7 _ i . 00 h'(^) = 0, the limiting distribution of maxfc 7™^ is of Gumbel-type. Note that this is true even when 
t < 1 which refers to heavy-tail distributions. Solving for 1 — F^ k {lx) = 7^, we get 

l K =P% 2a \ 2 logt^L. (11.29) 
Now, we apply the following theorem by Uzgoren. 

Theorem II.l (Uzgoren). Let x\ 1 . . . ,xk be a sequence of i.i.d. positive random variables with continuous and 
strictly positive pdf f x (x) for x > and cdf represented by Fx{x). Let hx(x) be the growth function. Then, if 
lim a; _ i . 00 g' (x) = 0, we have 

,.2 3 / -u+0(u 2 h' x (l K ))\ 

log { - logF*^ + h x (l K ) u)} = -u+ -h' x {l K ) + -(h x {l K )h' x {l K ) - 2h' x (l K )) +0( J . 

Proof: See [Uzg54, Equation 19] for proof. ■ 
The above theorem gives taylor series expansion of the limiting distribution for Gumbel-type distributions. In 



particular, setting Ik = /3 r a X 2 log* and u = loglogif, we have Ii(Ik) = 0[ — — ), H'Qk) 

V log * K J 



^( log if )' W'Qk) = — 2^ — ), and so on. In particular, we have 



log^ log K 1 
log K 



Pr(max^ k <l K + h(l K )loglogK) = e - e ~ IoeIoeK+ <H ».« , (n3Q) 

= 1-0(W), (11.31) 



\ogK 



where we have used the fact that e x = 1 + 0{x) for small x. Similarly, 



log log K + O ( l0! fj 6 K ) 

Pr[maxrt k <l K -h{l K )lo g logK) = e' e V g ' (11.32) 



A- 



(11.33) 

= °Q- (IL34) 

Subtracting (11.34) from (11.31), we get 

Pr(l K - loglogiC < max 7^ < l K + log log A') > 1 - . (11.35) 
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Note that the above equation is the same as (1.50). Therefore, following (I.51)-(1.61), we get 

and 



C* = BN O flog log 2 /* Kr2 ° 



C* = BNf™(r,B, ^^(loglog 2 / 4 ^) 



Further, if 



log 2/t KN 



N 



> 1, then C* =0[ log 



N 



(11.36) 
(11.37) 



C. LogNonnal 

In this case, \u" k \ is distributed according to LogNormal-(a, w) distribution. Hence, |^™ fc | 2 is distributed according 
to LogNormal-(2a, 4u>) distribution. The cumulative distribution function of 7™ fc , i.e., F^ k (j) (when 7 > 0) is 

tfyfcfr) = [p(WT,k\ 2 <^)fG i , k {9) d 9 Cn.38) 



= 1 - - 
2 



fa- 



erie 



log J- 2a 



fG Z:k {g) d 9, 



(11.39) 



//3 2 (p+d)- 2 ° 

where erfc[ ] is the complementary error function. Using the asymptotic expansion of erfc[-], F 7 « fc (7) can be 
approximated [AS70, Eq. 7.1.23] in the large 7-regime as: 



1 /•' 3 r ° 

F 7? fe (7) « 

2 7 /3 2 (p+d)- 2 ° (--^ ) 



/ /C?i,*(g) / W I -2a \ Z^( _1 ) 



(2m- 1)!! 



r d 5 , (11.40) 



where (2m — 1)!! = 1 x 3 x 5 x . . . x (2m — 1). We ignore the terms m — 1.2,... as the dominant term for large 
7 corresponds to m = 0. Therefore, we have 



l-i/ — 



o2 — 2q 

2w Z"' 3 r ° 



/3 2 ( P +d)- 2 ° log 2a 



2 
9 

— 3a N 2 



-fGi, k (g)dg 



(11.41) 



1 - 



2w rg e 



7r p 2 log l 2n - 2a 

P 2 r 



1 



3 



22 a/3 2 p 2 V/3 2 



log ^ -2c 



log 2 - 2a 



2w /" {p-d)*<* e 

(p + d) 2 °> 



log 2 - 2a 



-ds( 5 ). 



Now, we claim that 



lim (l-i^(7))(log 7 -log(/3 2 r - 2Q ) -2a)e 



rl I2w 



(11.42) 



(11.43) 
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This is because the contribution of the two integrals in (11.42) towards the RHS of (11.43) is zero. The contribution 
of first integral, when 7 is large, is 

( ' ) 



(log 



2a e 



< lo 



2^-Za 



-2a-2 



a(3 2 p 2 



log £- 2a 



-dp 



< 



< 



-2a-2 



a/3 2 p 2 



log 2 - 2a 



^3 



-dg 



' 



< 



aj3 2 p 2 

r -2a-2 

a(3 2 p 2 

-2a-2 
10 

afP-p 2 

-2a-2 



^^^^ 







( P -d) 
"5' 



S57 1°8" 



4(i ) 



0V 
3 



(p-d) 



/SV** 



dg 



a ? 2 P 2 ^(log^=^-4a 



1 - 



p — d 



— > 0, as 7 — > 00. 



(11.44) 
(11.45) 
(11.46) 
(11.47) 

(11.48) 

(11.49) 
(11.50) 



where in (11.44), we take an upper bound by taking the term (-p-) 1 out of the integral, and in (11.48), we 



put g = /3 r Q in the exponent of ( — 9 _ ia ) since g < P 2 r za . The second integral has an exponent term that 

V P' r 



goes to zero faster than e 



log ■ 



0, making its contribution zero. Note that only the first two term in 
(11.42) contribute to the RHS in (11.43). Similar to the above analysis, it is easy to show that 

:-2o\ 2 



1™ fi?Si)i e 







7—^00 



P 2 \ /, 8wit 



Using the above equation and (11.43), we have 



/i(7j = — j—^ — ps — for large 7, and 



/7",(7) log 7 



lim h' (j) 

7—^00 



Therefore, the limiting distribution of max^ 7™ fc belongs to the Gumbel-type. Solving for Ik, we have 



(11.51) 

(11.52) 
(11.53) 

(11.54) 
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KIk) = O(i^fe), h'(l K ) = O(^), //'(/k) = o( I - 1 i i7 -), and so on. Using Theorem H.l for u = 
log log A", we have 

Pr ( max 7 " fc < ^ + log log A) = e~ e (11.55) 



where we have used the fact that e x = 1 + 0{x) for small x. Similarly, 

loglogg + of 1 "^ 1 ^ ) 
p \ \/log K / 



Pr( max7 l r ; fc </ x -/i(/ K )loglogA') = e~ e (11.57) 



= P -(i+o(™)) log x 



e -\^\^nr)r^~ (11.58) 
- 1 
.A 



0(4)- CH.59) 



Combining (11.56) and (11.59), we get 

, e V8w log K e V&v) logK . /In 

Pr ( Za- - c , log log A < max 7 £ fc < Z K + c— =_ log log A J > 1 — Of — - , (11.60) 
V v log A fe ' v log A / V log A / 

where c is a constant. Now, following a similar analysis as in (I.51)-(1.61), we get 

max7jfc = QQk) w.h.p., (11.61) 



t/log^-). 



Further, if — ^ — > 1, then C* = O I log - 

III. Proof of Theorem 4 and Corollary 3 for regular extended networks 

The proof of Theorem 4 and Corollary 3 under the setup of regular extended networks is similar to the proof of 
Theorem 2 and Corollary 1 with the change of replacing p 2 by BR 2 . Here, we use Lemma 1.2, and Lemma 1.3 with 
p 2 = BR 2 . Lemma LI, however, is replaced by Lemma III. 1 (given below) wherein we show that the interference 
term in SNR expression can be bounded by a constant for extended networks. 

We now give the proof details. The counterpart of Lemma LI, i.e., Lemma III. 1, for regular extended networks 
is as follows. 

Lemma III.l. In a regular extended-network and Rayleigh fading channel, i.e., v™ k ~ £/V"(0, 1), the lower bound 
on achievable sum-rate of the system is 

C* > .^ N (r,7V)^E{log(l + P C onmax 7 ™ fc )}, (IH.l) 

i,n 

where r > is a fixed number, f^ N (r^ N) = ^^^ rcr ^ Co , [i and a are the mean and standard-deviation of |^ n /J 2 , 

32 2-2q 
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Proof: We have, from Theorem 1, 

c*>Ve{ maXfclQ g( 1 + p -^) 1 (m . 2) 

where the expectation is over user-location set (x,y) = {{xk,Uk) V /c}, and the fading random-variable set v = 
{ff k V i,k,n}. We will first consider the expectation w.r.t. v for a given (x,y). Then, ijj j. is known for all 
{j,k). We now apply one-sided variant of Chebyshev's inequality (also called Cantelli's inequality) to the term 
in the denominator, Ylj^i ^Jk^l^fkl 2 - Let the mean and standard deviation of \v™ k \ 2 be /i and a. Note that u 
is independent of (x, y) by assumption. Therefore, the mean and standard deviation of Y^jM . H-ik"^^ 2 are 
RJk a anc l a \Jj2j^i res P ec tively. Hence, applying Cantelli's inequality, we have 

1 



=► ^(Ek*i 2 >(^^)E^ 20 ) < 



1+1 
1 



1 + r 2 

.2 



where r > is a fixed number. 

Now, we break the expectation in (III. 2) into two parts — one with Y^j^ti Wjkl 2 > (/■* + rcr ) Sj^i ^7^" anc ^ 
other with l^fcl 2 — ra ^^39 L i ^Jk" ■ We tnen ignore the first part to obtain another lower bound. 

Therefore, we now have 

c , > r 2 " " E {max, log (1 + Peon 7l? fc )} 

Applying coordinate geometry in order to upper bound ^Jk*' me coordinate of any base-station can be 

written in the form of m(2R, 0)+n(R, i?v / 3), where m, n are integers, and 2R is the distance between neighbouring 
base-stations as shown in Fig. III. 1 . 
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We can write 

T, R i,t a=r 2aB <ro+ E R jl a - ( IIL5 ) 

where B< ro are the number of base-stations within ro distance of user k with corrdinates (xk,yk)- Now, 

-2a 

( \J (2mR + nR- x ) 2 + (nRVs - y ) 2 ] 

Rj,k>ro 



E = £ ( V(2™i? + ni?-xo) 2 + (ni?V3-yo) 2 



< 



d(x,i/)>r } / ^/ ( 2a; # + yiJ - xo) 2 + {yRs/l - y ) 2 

V v ' 

1 dy'dx 



(111.6) 



VT^m I I r TW drde (IIL7) 
2V3i? 2 Jr=r J 0=0 [r 1 ) 



{x«+v«>^} 2V3i? 2 (x' 2 + 2/' 2 ) c 

/■OO />27T 



2-2a 
7T r n 



V3i? 2 2a- 2' 

where in (III. 6), we substituted x' = 2xR+yR—Xo and y' = yRy/3 — yo, and in (III. 7), we substituted 

r 2 4r 2 

and y' = rsin#. Moreover, B< ro rj < Therefore, from (III. 5), we have 



(III.8) 



^,2-20 



E^fe a<i V- 4 +^7^ 7 • ( IIL9 ) 



Combining the above equation with (III.4), we get 



c . > EE fT^ m ax fc log(l + Pcon7i: fc )l 



P fl 2 r 2_2Q 



where c = „2° 4 + -=-^ — - is a constant. From (III. 10), we have 

n \ V3(2a— 2) J 

C* > /,f(r, iV) £e{ log (1 + max 7 " fc )}, (HLll) 

where /^(r, iV) = g^ . - 

Note that, for Rayleigh fading channels, 

|^™ fc | 2 - exponential (1) V i,k,n. (111.12) 

Therefore, \i = a = 1. The rest of the steps in the proof remain same as those in the proof of Theorem 2 and 
Corollary 1 in Section I with p 2 = BR 2 . In particular, applying Lemma 1.2, and Lemma 1.3 for extended networks 

if r 2 

yield Ik = /3 2 rQ a log -g^z- Therefore, using [SH05, Theorem 1], we get 

Pr{l K - log log ^ < max 7^ < Z* + log log ^ } > l-0(£^). (HL13) 

^g 7i? 

Note that the only change in the above equation from (1.50) is that K is replaced by K/B. Following the analysis 
in (I.51)-(1.61) (after replacing K by K/B), we get the final results in Theorem 4 and Corollary 3. 
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IV. Proof of Theorem 5 and Corollary 4 for regular extended networks 

For proof of this theorem and corollary, replace p 2 in the proof of Theorem 3 and Corollary 2 in Section II 
by BR 2 and use Lemma III. 1 (instead of Lemma 1.1) to obtain the lower and upper bounds for regular extended 
networks. 
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